arXiv:1506.03698vl [cond-mat.mtrl-sci] llJun2015 


Giant Magneto-Seebeck Effect in Spin Valves 


X. M. Zhang, C. H. WanQz. H. Yuan, H. Wu, Q. T. Zhang, X. Zhang, B. S. Tao, C. Fang, and X. F. Harfl 

Beijing National Laboratory for Condensed Matter Physics, 

Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
(Dated: June 12, 2015) 

Giant magneto-Seebeck (GMS) effect was observed in Co/Cu/Co and NiFe/Cu/Co spin valves. Their Seebeck 
coefficients in parallel state was larger than that in antiparallel state, and GMS ratio defined as (Sap-Sp)/Sp 
could reach -9% in our case. The GMS originated not only from trivial giant magnetoresistance but also from 
spin current generated due to spin-polarized thermoelectric conductivity in ferromagnetic materials and subse¬ 
quent modulation of the spin current by different spin configurations in spin valves. Simple Mott two-channel 
model reproduced a -11% GMS for the Co/Cu/Co spin valves, qualitatively consistent with our observations. 

The GMS effect could be applied simultaneously sensing temperature gradient and magnetic field and also 
be possibly applied to determine spin polarization of thermoelectric conductivity and Seebeck coefficient in 
ferromagnetic thin films. 


Spin caloritronics, focusing on interplay between thermo¬ 
electric transport and spin configurations of ferromagnetic 
materials or related hybrid structures, has extensively re¬ 
shaped landscape of recent spintronics research, since dis¬ 
covery of spin Seebeck effect)IHg], spin-dependent Seebeck 
effect1213,11 and spin-dependent tunneling Seebeck effect(9f- 
HHl . The former two effects have shed light on effective meth¬ 
ods of generating a steady pure spin current into nonmag¬ 
netic materials through wasted heat which is being abundantly 
produced and wasted in current microelectronics. Thereafter 
marriage of thermally generated spin current and varied spin 
configurations in magnetic tunnel junctions has bred spin- 
dependent tunneling Seebeck effect. However, another ver¬ 
sion of Seebeck effect in broadly applied spintronic structure, 
namely, giant magnetoresistance (GMR) spin valves, which 
deals with spin-configuration-dependent thermoelectric trans¬ 
port in diffusive manner, has still not been systematically re¬ 
searched. Gravier|J_2] and Sakurai|13ll et al. reported a field- 
dependent thermopower in Fe/Cr and Co/Cu multilayers and 
they discussed little on physics hidden behind it. Very re¬ 
cently, Jainj ll et al. reported a tiny field-dependent Seebeck 
effect in CoFe/Cu/CoFe/IrMn/Ru spin valves. They further 
attributed its field dependence to the trivial origin that differ¬ 
ent resistances in parallel (P) state and antiparallel (AP) state 
resulted in different Seebeck coefficients accordingly. Fur¬ 
thermore, HulfH et al. reported a small field-dependent See¬ 
beck coefficient in nonlocal lateral spin valve which showed 
even different sign of (Sap-Sp)ISp with other literatures [12, 
13]. Magneto-Seebeckratios [GMS=(S'ap-*S'p)/>S'p] reported 
in these researches not only diverged in sign but also in abso¬ 
lute values ranging from 0.1% to 10%. In this letter, we will 
demonstrate that Seebeck effect of spin valves not only ex¬ 
hibit specific field dependence as that of GMR effect but also 
its GMS ratio could be giant and even reach a higher value 
than GMR ratio. Furthermore, the GMS ratio has different 
sign with GMR in our case, which could no longer be ex¬ 
plained as a parasitic effect of the GMR effect. We will also 
show physical origin behind this phenomenon in this letter af¬ 
terwards and give a probable explanation on sign change in 
different material systems. 


As shown in Fig. 1(a), spin valves with structure of 
Si0 2 //Ni 8 iFei g (3)/Cu(f)/Co(4)/Ir 22 Mn 7 8(12)/Pt(3) (nominal 
thickness in nanometer) were grown by ULVAC magnetron 
sputter system with a base pressure of l.Ox ICG 6 Pa at room 
temperature. During deposition, a field of 200 Oe was used 
to induce an easy axis (EA) as well as exchange bias (EB) 
along x in the spin valves. Different samples with f=2, 4, 6, 
8 correspond to Sample 1-4. All films were then patterned 
into 20 fj,mx210 fim bars as shown as the dark green bar in 
Figure 1(b) by ultraviolet lithography and subsequent argon 
ion milling. Then a bar made of Cu(20)/Au(20) between Pad 
C and Pads D as well as Pad A was deposited to transport 
heating current and detect Seebeck voltages, respectively. As 
measuring magneto-Seebeck effect, we applied a heating cur¬ 
rent J=/osin(u;i) by Keithley 6221 between Pad C and Pad D 
along y which produced a thermal gradient along x (V :/: 7’). In 
this direction, no net current flowed between Pad A and Pad 
B since the bar made of spin valve (SV) were isolated from 
the heating bar by 7 //m. Pad A and Pad B were just used to 
measure Seebeck voltages along the SV. Seebeck voltage was 
first amplified by a low-noise preamplifier (SR560, Standford) 
with gain of 20000. Then SR830 only measured locked-in 
voltage signals 14s=14socos(2cuf) to improve signal-to-noise 
ratio and also to rule out voltage origins such as Peltier ef¬ 
fect and other higher order signals. Magnetic field was pro¬ 
vided by Physical Properties Measurement System (PPMS- 
9T, Quantum Design). Environment temperature during mea¬ 
surement was kept constantly as 300 K by PPMS. We have 
verified that GMS was indeed independent on oj though the 
magnitude of Vso did (SI A). In the following, we will only 
present results measured at w=731.7 Hz. M-H hysteresis 
loop was measured by vibrating sample magnetometer (EV9, 
MicroSense). 

Fig. 2(b) shows mainly the dependence of Hgo on heat¬ 
ing current Iq for Sample 3 (t= 6). Fig. 2(a) shows the 
field dependence of two-terminal resistance between Pad A 
and Pad B as applied B along x. Exchange bias of the 
Co pinned layer could be clearly resolved with a bias field 
of about 200 Oe. The hysteresis loop around zero fields 
originated from the NiFe free layer. Fig. 2(a) also shows 
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FIG. 1. (Color online) (a) structures of spin valves and (b) pattern of 
devices and second harmonic measurement setups. 


a natural result that resistance of Sample 3 in antiparal¬ 
lel state (AP) was larger than that in parallel state (P), as 
normal GMR sandwiches. Interestingly, very similar to 
the field dependence of resistance, a field-dependent See- 
beck voltage was also observed. However, Seebeck volt¬ 
age in AP state was remarkably smaller than that in P 
state. Here we have defined GMR(B)=[p(B)-p(0)]/ p(500 
Oe) and GM S(B)=[S(B)-S(0)]/S(500 Oe)=[V S o(B,I 0 )- 
Hso(0,/o)]/^So(500 Oe,/ 0 ) with p(B) and S(B ) being resis¬ 
tivity and Seebeck coefficient at B. For Sample 3, GMR (500 
Oe)=1.5% while GMS {500 Oe)=-9%. The value of GMS 
was not only much higher than that of GMR, but also much 
higher than other thermoelectric signals such as anisotropy 
Seebeck effect as shown in the following results. Therefore 
this magneto-Seebeck effect was named as giant magneto- 
Seebeck effect. Fig. 2(c) shows that Seebeck voltage in P state 
linearly depended on heating power ( P ) as expected while the 
value of GMS (500 Oe) was nearly independent on the P. It 
was also significant that exchange bias field of the Co layer 
was gradually reduced due to the increase in average temper¬ 
ature of the sample. Higher temperature usually results in 
higher switching probability in the same held due to higher 
thermal perturbation as reported in Ref. [ Thi F7 ]. 


In order to rule out other possibilities of the field-dependent 
voltage contribution from anisotropy Seebeck effect under 
V a T or anomalous Nernst effect under V.T or planar Nernst 
effect under V y 7’, we have performed some controlled ex¬ 
periments. Fig. 3(a) shows M-H hysteresis loop of Sam¬ 
ple 1(1=2) and Sample 2(1=4) with B along the EA and 
EB direction. Different from Sample 2, pinned layer and 
free layer in Sample 1 coupled with each others. There¬ 
fore only P state in Sample 1 was observed, which agreed 
with magnetotransport measurement. We only observed a tiny 
anisotropy magnetoresistance as well as a planar Hall volt¬ 
age in Sample 1 as varying B along x, without any signals 
from GMR effect (Fig. 3(b) and (c)). Meanwhile, as we 
applied / () =80 mA along the heating bar and B along the 
hard axis ( y ), only a negligible small field-dependent volt¬ 
age was observed around zero fields (Fig. 3(b) and (c)). 
The magnitude of the field-dependent voltage was only about 
15 nV, close to the noise level of 10 nV but much smaller 
than that from GMS effect. The latter reached over 100 nV 
as / (l =80 mA in Sample 3. This tiny field-dependent See- 
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FIG. 2. (Color online) (a) two-terminal resistance between Pad A and 
Pad B in Sample 3, (b) field-dependent Seebeck voltages measured 
at elevated heating current Jo and (c) the dependence of Seebeck 
voltages in parallel state on heating power P=I%R\ l where f?h=80 
was resistance of heating bar. 
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FIG. 3. (Color online) (a) M-H curves of Sample 1 and Sample 2. 
The latter shows separated switching of free layer and pinned layer 
in Sample 2 (t=4) while the former shows parallelly coupled free 
layer and pinned layer in Sample 1 (f=2). Anisotropic magnetore¬ 
sistance (black and empty squares) and thermoelectric voltages (blue 
and empty circles) of Sample 1 as (b) B along easy axis and (c) B 
along hard axis. 


beck effect in Sample 1 was induced by anisotropy magneto- 
Seebeck effect. We also measured field-dependence of See¬ 
beck voltages of Co(20)/Pt(3) and NigiFeig(20)/Pt(3) and 
also found no observable field-dependent voltages (SI B). 
These controlled measurements compelled us to conclude 
that the field-dependent voltages in Fig. 2 and Fig. 4 re¬ 
sulted from different Seebeck voltages between AP state and 
P state, instead of from planar Nernst effect or anisotropy 
magneto-Seebeck or anomalous Nernst effect. GMS and 
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FIG. 4. (Color online) (field dependence of GMS and GMR4 (mea¬ 
sured in four-terminal geometry) in (a) Co/Cu/Co/IrMn/Pt and (b) 
NiFe/Cu/Co/IrMn/Pt spin valves 


GMR effects were also detected in another spin valve 
Si0 2 //Co(4)/Cu(4)/Co(4)/Ir 22 Mn 78 (12)/Pt(3) (Fig. 4). GMR 
and GMS ratios of the Co/Cu/Co spin valve were, respec¬ 
tively, 4.6% and -6.0%, both higher than 2.6% and -4.5% of 
Sample 2. Though different in values, these spin valves shared 
some characteristics in magneto-electrotransport properties 
and magneto-thermotransport properties. Firstly, their GMR 
ratios were all positive as usual while their GMS ratios were 
always negative. Results measured at low temperatures also 
showed the same trends (SI C). This universal phenomenon 
was also observed in Sample 4 (not shown here). The negative 
GMS ratios meant that |S A p| < |Sp|. Secondly, their GMR 
ratios were all smaller than their GMS ratios in absolute val¬ 
ues. In the following discussion, we will introduce physical 
origins beneath this universality based on a simple linear re¬ 
sponse theory in which temperature gradient and electric field 
independently influence charge and spin transport inside SV 
and in which they do not affect transport coefficients such as 
conductivity and thermoelectric conductivity of the system. 
Mott two-spin-channel model[ 18-2(|| which assumes much 


longer spin relaxation length than mean-free path in both fer¬ 
romagnetic and nonmagnetic layers has been broadly applied 
in estimating GMR values. Here we generalized the model 
into thermotransport region and tried to evaluate the value of 
GMS. For spin relaxation length in copper is several hundreds 
of nanometers 112114231. much higher than thickness ( t ) of cop¬ 
per in our spin valves, we further supposed there was nearly no 
spin relaxation event occurring in the copper. We also suppose 
the same thickness of the two FM layers as occurred in our 
Co/Cu/Co spin valves. Besides, we theoretically dealt with the 
case of current-perpendicular-plane (CPP) geometry instead 
of current-in-plane (CIP) geometry to grasp main physics 
though practical measurement was implemented in current- 


and temperature-gradient-in-plane geometry. In P state, spin- 
polarized current in the first FM layer would equal that in the 
second FM layer for each spin channel since few spins flip 
inside Cu as assumed. J^=al A VT/ A =a| A VT 2 t/ WJ /J ' 
to fulfill continuity condition of current density in both spin 
channels. Here jJ A , and VXf A were spin-polarized 
current, thermoelectric conductivity and temperature-gradient 
in zth FM layer for t/| spin channels, respectively. Es¬ 
pecially, the concept of different temperatures in different 
spin channels has also been supposed and even experimen¬ 
tally testified in Ref. [250- Total current J=j\+ 4=jJ+J 2 
and average temperature gradient 2VT=VT/ A +VT 2 A , ne¬ 
glecting temperature gradient dropped on copper. Thus J = 
+ ct 2 ) + a\a\/{a\ + a 2 )]2VT = 2apVT. 
The ap is thermoelectric conductivity of spin valves in 
P state. In this case, pp={p\+pl){p\+p^)l{p[+p\+p[+p^). 
Therefore Sp = Ep/VT = ppJ/VT = 2ppap. 
In AP state, jJ /l = 4 A VT 1 TA = a At VT 2 J ' /t = 
J 2 A also to satisfy current density continuity requirement. 

J = 4 + 4 = 4 + 4 and 2VT=VT 1 1 ' A +VT 2 J7t . 

J=[44/(4+4)+at4/(at+4)]2VT=2a A pVT. Here, 
correspondingly, the a A p is thermoelectric conductivity in 
AP state. And pAP=(Pi+P 2 )(Pi+P 2 )/(pi+P 2 +Pi+/4)- Thus 
S\p=Ep L pl'S/T=p\pJI'S/T=2p\paxp- Ratio of GMS could 
be expressed via EquationQ] 

GMS + 1 = (P | + pi f p \ + P Y a \ + 1 )( 1 + ( 1 ) 

(p! + p\){Pi + piM + i + a\) 

It shows that GMS ratio is dependent both on spin polariza¬ 
tion of resistivity and thermoelectric conductivity in two FM 
layers. For the Co(4)/Cu(4)/Co(4) spin valve, GMR and GMS 
could both be simplified as Equation[2]and Equation^ Equa- 
tion [2] is also derived in Ref. (26] and it is still applicable for 
CIP geometry. We thus expect Equation [2 is also instructive 
enough to explain our thermotransport data in CIP geometry. 

GMR =Y^p2 (2) 


GMS = 


p2 _ p2 

r (T 1 ft 

i -Pi 


(3) 


Here P CT =(<T t -er i )/(er t +cr i ) and P Q =(a r a 4 .)/(a t +a 4 .), cr t /4 

being conductivity of spin-up/down channel of FM layer. 
a t/I =(T t/4 5 1'/4- Thus P a =(P s +P a )/(\+P s P a ) and P s =(S r 
Si)/(S^+Si) where is Seebeck coefficient for spin- 

up/down electrons in FM. The above derivation could also 
be interpreted in terms of a generalized equivalent circuit as 
shown in Fig. 5. V and VT independently determine drift¬ 
ing and diffusion current density inside SV, respectively. To 
the role of p in ordinary electrotransport theory, reciprocal of 
thermoelectric conductivity ( p/S ) plays a similar role in ther¬ 
moelectric transport theory. The equivalent circuit in Fig. 5 
could reproduce the results of Equation |T][2 
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FIG. 5. (Color online) equivalent circuits driven by W in (a) AP 
state and (b) P state and also driven by VT in (c) AP state and (d) P 
state, respectively. 


Equation [2] and [3] shows that always GMR >0 while 
GMS <0 as 
P n ~ 0.5 


P a | < \P a \. For Co, P(j = 0.4 [22] while 
2%, thus ideal GMS ~-ll%<0 according 
to Equation[3] semi-quantitatively consistent with our experi¬ 
ments. |Sp| is also larger than | S'ap | in the research of Walter 
@1. If P s =0, P a =P a . Thus GMS= 0. The nonzero GMS ef¬ 
fect in our research indicated that Seebeck coefficient in cobalt 
was indeed spin-dependent as reported in Ref 1231 24]. If 
P a = 0, GMR=GMS. The observation that GMR and GMS 
had different signs in our spin valves also indicated that a in 
Co was spin-polarized. If P a = 0, GMR =0 while GMS=-P%, 
meaning GMS effect could still survive even without GMR. 
This specialty is very different from that derived in Ref. [12]. 
The above derivation shows it was (1) different thermoelec¬ 
tric conductivities in different spin channels and (2) continuity 
condition of current density in both spin channels and also in 
both AP and P state that led to a giant magneto-Seebeck effect 
in our spin valves. 

In summary, giant magneto-Seebeck and giant magnetore¬ 
sistance effect were unambiguously observed in Co/Cu/Co 
and NiFe/Cu/Co spin valves. The former effect led to a neg¬ 
ative GMS value which could reach -9% in experiment and 
be estimated as about -11% in a simplified two-spin-channel 
model. The GMS value was much higher than anisotropy See¬ 
beck effect and other possible parasitic signals. This GMS 
effect was attributed to nonzero spin polarization of Seebeck 
coefficient of FM layers and modulation of spin-polarized cur¬ 
rent in different spin configurations. This research not only 
hinted a prospective application of spin valves in high sen¬ 
sitive temperature and magnetic field sensors but also in T 
driven field-sensors. GMS also provided a probable manner 
to estimate spin polarization of Seebeck coefficient of a cer¬ 
tain FM material which was very hard to measure in current 
stage. 
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FIG. 1. (Color online) Field-dependence of Seebeck voltages mea¬ 
sured in different frequencies. 
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FIG. 2. (Color online) Field-dependence of GMS ratios measured in 
different frequencies. 



A. Frequency dependence of Seebeck effect of 
Co/Cu/Co/IrMn/Pt spin valves 

We have measured frequency dependence of Seebeck volt¬ 
ages and GMS ratios. As shown in Figure SI, except results 
measured in 3.1 Hz, the Seebeck voltage in parallel state were 
gradually decreased with increasing frequency in our AC mea¬ 
surement, which might be probably induced by a lower heat¬ 
ing efficiency in higher frequencies of heating current. How¬ 
ever, GMS ratios were insensitive to the frequencies as shown 
in Figure S2, hinting this ratio was an intrinsic property of our 
structure or materials. No relevance between the GMS ratio 
and the frequency as expected also indicated effectiveness of 
our measurement geometry and 2nd harmonic technique. 

B. Measurement of Seebeck effect of thick NiFe and Co 
films 

To further confirm that giant magneto-Seebeck (GMS) ef¬ 
fect was not due to anisotropic magneto-Seebeck effect of 
NiFe and Co films, we also measured Seebeck effect in NiFe 
and Co thick films. They were both independent on field 
within our measurement accuracy as shown in Figure S3. Be¬ 
sides, the measured ratio of SmFe/Sco was close to 1, also 
close to those literature values of bulk NiFe and Co JT|]. 

C. T- and if-dependence of GMR and GMS ratios of 
Co/Cu/Co spin valves 

Figure S4(a) and Figure S5(a) shows that GMR ratios in¬ 
creased linearly with decreasing temperatures and the GMR 
ratio increased from 5.5% at 300 K to 9% at 50 K while ab¬ 
solute value of GMS ratio decreased with decreasing temper¬ 
atures, from 6% at 300 K to 0% at 50 K as shown in Figure 
S4(b) and Figure S5(b). Same as the results at room temper- 
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FIG. 3. (Color online) Field dependence of Seebeck coefficients of 
thick NiFe and Co films. 
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ature, GMR > 0 while GMS < 0 at low temperatures as 
well. dGMR/dT=-l.5%/l00 K, as occurred in other GMR 
spin valves 0. Regarding saturated GMS ratios at high tem¬ 
peratures, we think it was related with energy-dependence 
of GMR ratios as following discussion. According to Mott 
relation jsj], S = So(k B T)(dliV7/dE)E F , So=Tt 2 k B /3e=283 
fiV/ K, where S, k B , T, a, E and E p is Seebeck coef¬ 
ficient, Boltzmann constant, Kelvin temperature, conductiv¬ 
ity, carrier energy and Fermi energy, respectively. GMS = 
(Sap — S P )/Sp. 

-GMS • S P = —GMS ■ S 0 (k B T)(d\na P /dE) EF 

=S P - S A p=So(fc B T)(dlncrp/dF;)-So(fcBr)(dlncrAp/dF;) 

=So(k B T)dln(ap/aAp)/dE=S 0 (k B T)d\n(l+GMR)/dE. 
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FIG. 4. (Color online)Field-dependence of GMR and GMS ratios 
under different temperatures 


In our case, GMR -C 1, therefore, dGMR/dE « 
—GMS xS P /(ksTSo), which indicated that GMS ratio had 
close relation with energy derivative of GMR ratio. S p /T in 
our experiment was nearly constant, about -0.015 V/K 2 , dur¬ 
ing 150 K-300 K. Thus dGM/f/d/(=-3.7%/eV, higher energy 
of carriers resulting in lower GMR ratios, as reported in Ref. 

II- 


* Electronic mail: wancaihua@iphy.ac.cn 
' Electronic mail: xfhan@iphy.ac.cn 

[1] F. K. Dejene, J. Flipse, and B. J. van Wees, Phys. Rev. B, 86, 
024436 (2012). 

[2] J. Garcia-Torres, E. Vails, E. Gmez, J. Magn. Magn. Mater. 322, 
3186 (2010). 

[3] M. Cutler, N. F. Mott, Phys. Rev. 181, 1336 (1969). 

[4] Y. Sato, K. Hoshino, S. Okamura, K. Kato, H. Hoshiya, IEEE T. 
Magn. 46, 1610 (2010). 



FIG. 5. (Color online)Temperature dependence of maximum GMR 
and GMS of Co/Cu/Co/IrMn/Pt spin valves. 


































